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ABSTRACT 


A Comparison of Experimental and Theoretical Results for Leakage, 
Pressure Distribution, and Rotordynamic Coefficients for Annular Gas 

Seals. (December 1984) 

Colby Oran Nicks, B.S., Virginia Polytechnic Institute and State 

University 

Chairman of Advisory Committee: Dr. Dara Childs 

This thesis concerns a study of annular gas seals which is 
currently in progress at Texas A&M University. A brief discussion of 
the importance of seal behavior in rotordynamics is presented, as is a 
review of current annular seal theory. An outline of Nelson's 
analytical-computational method for determining rotordynamic 
coefficients for this type of compressible-flow seal is included. 
Various means for the experimental identification of the dynamic 
coefficients are outlined, and the method employed at the TAMU test 
facility is explained. The TAMU test apparatus is described, and the 
test procedures are discussed. Experimental results, including leakage, 
entrance-loss coefficients, pressure distributions, and rotordynamic 
coefficients for a smooth and a honeycomb constant-clearance seal are 
presented and compared to theoretical results from Nelson's analysis. 
The results for both seals show little sensitivity to the running speed 
over the test range. Agreement between test results and theory for 
leakage through the seal is satisfactory. Test results for direct 
stiffness show a greater sensitivity to fluid prerotation than 
predicted. Test results show that the deli berately-roughened surface of 
the honeycomb seal provides improved stability versus the smooth seal. 
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NOMENCLATURE 

A,B = Fourier coefficients for rotor motion 
C,c = direct and cross-coupled stiffness coefficients (FT/L) 
e 0 = displacement of seal rotor from centered position (L) 
K,k = direct and cross-coupled stiffness coefficients (F/L) 
k = entrance-loss coefficient 

M,m = direct and cross-coupled added-mass coefficients (M) 
m = fluid mass flow rate (M/T) 

ns ,ms = stator Hirs coefficients 

nr,mr = rotor Hirs coefficients 

p = fluid pressure (F/L 2 ) 

R = seal radius (L) 

R a = 2pUC/p = nominal axial Reynolds number 
U = mean fluid flow velocity (L/T) 

X ,Y = radial seal displacements (L) 

Y = ratio of specific heats for air 

e 0 = e 0 / C r = equilibrium eccentricity ratio 

p = fluid density (M/L 3 ) 

A = Fanning friction- factor 

t = fluid shear stress (F/L 2 ) 

w = shaft angular velocity (1/T) 

Q = shaft precessional velocity (1/T) 

u = fluid viscosity (FT/L 2 ) 
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INTRODUCTION 

With turbomachinery design trends tending toward increased speeds 
and loadings, lighter weight, and reduced clearances between rotating 
and stationary parts, considerable concern with instability and 
synchronous response has arisen. Synchronous response refers to 
vibration of the turbomachine rotor assembly at a frequency coincident 
with the rotational speed. Characteristically, the vibration amplitude 
increases to a maximum at each critical speed (coincidence of the 
running speed with a rotor's damped natural frequency), and then 
decreases to a relatively steady level. Operation of turbomachines at 
rotational speeds above any of the critical speeds requires the rotor to 
traverse them during start-up and shut-down. Therefore, in order to 
limit the peak synchronous vibration levels, the machine designer 
aspires to introduce damping into the rotor system. 

In contrast to synchronous vibration, "unstable" or "self-excited" 
motion is typically subsynchronous . This motion takes the form of 
whirling of the rotor shaft at a natural frequency less than the 
rotational speed. The exciting force for this whirling motion is a 
tangential force acting on the rotor due to some fluid or friction 
mechanism. This vibration often occurs with large amplitudes which 

sustain or grow as running speed increases. At best, this self-excited 
whirling prevents further speed increases; at worst, it results in 
damage to, or catastrophic failure of the equipment. One of the 
rotordynamic force mechanisms which plays a role in self-excited rotor 
motion and synchronous response is that of the forces developed by 
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annular seals. Until recently, most investigations of annular seals in 
turbomachinery have been concerned with reducing the leakage of the 
working fluid through the seal (i.e., improving the sealing effect). 
However, recent experiences have shown that forces developed by these 
seals can have considerable influence on the stability and synchronous 
response of rotating machinery. Black et al . [1-3] have demonstrated 
the critical effects that forces developed by neck-ring and interstage 
seals have on the rotordynamic behavior of pumps. Also, stability 
difficulties with the Space Shuttle Main Engine (SSME) high-pressure 
fuel turbopump [A] have prompted further research into these forces 
developed by liquid seals. 

Experiences have shown that various gas seal configurations can 
have similar influences on the rotordynamic behavior of turbomachinery . 
In the high-pressure oxygen turbopump of the SSME, for example, initial 
vibration problems were remedied by changing the turbine interstage seal 
from a stepped-labyrinth configuration to a convergent taper seal with a 
honeycomb stator and a smooth rotor[5]. A lack of experimental data to 
completely explain this and other gas seal behavior makes obvious the 
need for research in this area. 

The purpose of this report is twofold. It describes the test 
facility and initial test program developed to experimentally measure 
the fluid forces induced by annular gas seals, and it provides a 
comparison of theoretically predicted and experimentally obtained data 
for smooth and honeycomb seals. The leakage of the working fluid 
through the seal, the pressure gradient along the seal length, entrance 
pressure-loss data, and rotordynamic coefficients provide a basis for 
comparison. A short discussion on seal theory is included, and various 
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rotor dynamic coefficient identification schemes are described. The work 
presented herein is intended to add to the rapidly expanding database on 
seal forces, and to determine the validity of one theoretical analysis 
for predicting those forces. 
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ANNULAR SEAL ANALYSIS REVIEW 

As related to rotordynamics , seal analysis has the objective of 
determining the reaction forces acting on the rotor arising from shaft 
motion within the seal. Due to similarities between plain journal 
bearings and annular seals, seal analysis is generally based on 
governing equations which have previously been developed for bearings. 

Annular seals and plain bearings are geometrically similar, but 
seals typically have radial clearance-to-radius ratios on the order of 
0.005, versus C r /R ratios of 0.001 for bearings. Due to seal clearances 
and pressure differentials, fully-developed turbulent flow normally 
exists. Also, seals are nominally designed to operate in a centered 
position. Journal bearings, on the other hand, have operating 

eccentricities which vary with running speed and load. Therefore, most 
of the rotordynamic work for bearings has been done to determine dynamic 
coefficient versus eccentricity relationships. 


Two linearized 

seal 

models , 
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terms 

of dynamic 

coefficients , 

have 
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relationship. 
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where the dynamic coefficients {KXX, KYY , CXX , CYY , MXX , MYY} and { KXY , 
KYX, CXY , CYX, MXY, MYX} represent the "direct" and "cross-coupled" 


stiffness, damping. 


and added-mass 


terms, respectively. These 
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coefficients are functions of the equilibrium eccentricity ratio 
c 0 = e o / C r , where the eccentricity ratio e Q equals the displacement 
(e 0 ) of the rotor from the centered position divided by the nominal 
radial clearance (C r ). The term "cross-coupled" refers to the coupling 
effect exhibited by the off-diagonal terms; specifically, motion in one 
plane introduces reaction forces in an orthogonal one. These 
cross-coupled terms arise from the fluid's circumf erential velocity 
component, and show a strong dependency on both the magnitude and 
direction (with respect to rotor rotation) of the velocity. This 
circumferential velocity component may arise from the prerotation of the 
fluid as it enters the seal due to some rotating element upstream, or it 
may develop as the fluid passes through the seal, with rotor shear 
forces "dragging" the viscous fluid around its periphery. The 
cross-coupled stiffness term usually produces a destabilizing force 
component , and therefore is of considerable interest. The cross-coupled 
damping and added-mass terms are generally much less influential than 
the cross-coupled stiffness term with respect to stability. For no 
fluid rotation, these cross-coupled terms are zero. 


The second linearized seal model applies for small motions of the 
rotor about a centered position in the seal, as shown in Fig. 2 . This 
model can be expressed 



where the dynamic coefficient matrices are skew-symmetric. 

Theoretical work on annular seals has been done for both 
incompressible and compressible fluids. Black et al . [6] have developed 
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analytical "short-seal" solutions for incompressible seals, which 
account for circumferential fluid flow due to wall shear stresses but 
not pressure perturbations . The analysis employs a bulk-flow assumption 
and accounts for fluid prerotation as is enters the seal. Childs' [7] 
incompressible seal analysis provides "finite-length " solutions, in 
which both shear and pressure- induced flow are' included. Childs' 
utilizes Hirs' [8] turbulent bulk-flow model, and accounts for inlet 
swirl as well as perturbations in axial and circumferential Reynolds 
numbers due to clearance perturbations. 

Compressible flow in seals has been analyzed by Fleming [9, 10] and 
Nelson[11, 12]. Fleming presents a short seal solution for the leakage, 
direct stiffness, and direct damping coefficients for straight and 
tapered, smooth, annular gas seals, but does not include the 
cross-coupled damping terms. Nelson, whose analysis is used for 
comparison in this report, analyzes both smooth and surf ace-roughened 
annular seals in the straight and tapered configurations. An outline of 
Nelson's analysis is included in the section that follows. 
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NELSON'S ANALYSIS 


Nelson [11, 12] has developed an analysis which provides both 

static and dynamic results for annular gas seals. The static results 
include fluid leakage through the seal, pressure gradient along the seal 
axis, and the fluid axial and circumferential velocities through the 
seal. Dynamic data provided by the analysis consists of the 

rotordynamic coefficients (direct and cross-coupled stiffness and 
damping terms) for small rotor motion about a centered position 
(equation(2) ) . Nelson assumes that the added-mass terms are negligible 
for gas seals, and, hence, equation(2) is written 



Nelson utilizes a modified Hirs' [8] turbulent bulk-flow fluid 
model to develop governing axial and circumferential momentum equations, 
and his model is completed by the continuity and energy equations. 
Hirs' model defines the wall shear stress t w as 

t w - 1/2 pU m 2 no(2pU m H / u) mo = 1/2 pU m 2 noR a mo (4) 
where U m is the mean flow velocity relative to the surface „upon which 
the shear stress acts, and H is the local seal clearance. Hirs' 
formulation assumes that the surface roughness is the same on the rotor 
and stator. However, if the bulk-flow velocities relative to the rotor 
and stator are substituted in equation (4), the shear stresses at the 


rotor and stator are, respectively, 
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(5) 


T r = 1/2 pU r 2 nr(2pU r H / y) mr 
Tg = 1/2 pUg 2 ns (2pU s H / p) ms 
Hence, different surface roughnesses in the seal elements can be 
accounted for via the empirical coefficients mr, nr and ms, ns for the 
rotor and stator surfaces. These coefficients may be calculated from 
static-pressure-gradient test data, and are then provided as input 
parameters for Nelson's analysis. 

Assuming small motion of the rotor about a centered position, 
Nelson uses a perturbation analysis similar to that employed by Childs 
[7] to develop zeroth and first-order perturbation equations. The 
zeroth-order solution represents a zero-eccentricity flow condition, 
with rotor rotation but without precession. This solution is iterative 
and yields the mass-leakage flow rate, and the axial distribution of 
pressure, axial velocity, density, and circumferential velocity. 

An iterative solution scheme is employed, using initial guesses for 
the zeroth-order seal entrance Mach number and entrance pressur e-loss 
coefficient. The entrance-loss relationship is defined by 

Y/(Y-1) 

p o (0) = 1 / {1 + [ (Y-1 ) (k+1 )M o 2 (0)] / 2} (6) 


where p o (0) is the seal entrance/reservoir pressure ratio and M o (0) is 
the entrance Mach number. The entrance Mach number is iteratively 
adjusted, and the loss coefficient k is recalculated according to a 
curve fit by Deissler [13] 

k + 1 =V 5.3 / log 10 R a (7) 

which is plotted in Fig. 3- At axial Reynolds numbers above 200,000, k 
is equated to zero. The iterative solution procedure for M o (0) and k 
continues until either: 



ENTRANCE LOSS curve 
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1) the Mach number at the exit reaches unity and the exit pressure 
is greater than the sump pressure (choked flow), or 

2) the exit pressure equals the sump pressure and the exit Mach 
number is less than unity (unchoked flow). 

The pressure, density, and velocity distribution and their 
derivatives which are determined in the zeroth-order solution are used 
in defining coefficients of the first-order perturbation equations. 
These equations define the pressure, density, and axial and 
circumferential velocity perturbations due to rotor motion, and are 
transformed to sixteen ordinary differential equations in the axial 
coordinate z. The four physical boundary conditions required for 
the solution of these equations depend on the perturbation conditions 
that are specified at the seal entrance and exit. The inlet 
circumferential velocity perturbation is zero. Expansion of the 

entrance pressure-loss relationship of equation (6) yields a second 
boundary condition. For choked flow, the first-order perturbation in 
the exit Mach number is zero, while for unchoked flow, the first-order 
perturbation in the exit pressure is zero. 

Application of these boundary conditions and numerical integration 
of the ordinary differential equations provides the first-order 
solution. Integration of the first-order pressure solution along and 
around the seal periphery yields the direct and cross-coupled stiffness 
and damping coefficients, K, k, C, and c, respectively. 

The input parameters which can be varied in Nelson's analysis 
include: 

1 ) reservoir pressure and temperature. 


2) sump pressure. 



3) seal geometry (i.e. radius, length, clearances), 

4) rotor rotational speed and precession rate, 

5) entrance circumferential velocity of fluid, 

6) rotor and stator surface roughness (Hirs constants), 

7) empirical entrance-loss relationship, e.g., Deissler's, 
equation ( 1 6 ) , and 

8) fluid viscosity, gas constant, and ratio of specific heats. 

It is apparent that a large amount of theoretical data can be 
generated to determine the influence that these various parameters have 
on the fluid forces in annular gas seals. However, there is a lack of 
experimental data with which to compare the results of Nelson's 
analysis. Currently, test results due to Wachter and Benckert [14] exist 
for labyrinth seals, a special class of non-contacting seals which have 
stepped surfaces or "teeth" on the rotor, stator, or both. Experimental 
results for smooth and/or surface-roughened gas seals are limited to 
data for honeycomb seals also published by Wachter and Benckert. Hence, 
the need for a test apparatus which can be used to study the effects of 
the same variables provided for in Nelson's analysis is obvious. The 
experimental data generated by such an apparatus would be valuable for 
comparison to both Nelson's theories and others which may be developed 


in the future. 



TEST CONCEPTS 


A number of test programs have been implemented to measure the 
stabilizing and destabilizing fluid forces which are developed by 
turbomachinery elements. Some are concerned mainly with the study of 
seal forces, while others examine the forces developed by centrifugal 
pump impellers. In each case, reaction force and relative motion 
measurements are used for rotordynamic coefficient identification. Four 
general approaches have been employed, and will be reviewed here. 

Wachter and Benckert [14] employ a static displacement method for 
determining stiffness coefficients. In this method, as shown in Fig. 4, 
the rotor is displaced statically to some measured eccentric position 
while a pressure differential forces the working fluid past the seal. 
By measuring the reaction force components which are parallel and 
perpendicular to the static displacement vector, the direct and 
cross-coupled stiffnesses can be determined. Referring to equation (2) 
for small rotor motion about a centered position, a static rotor 
displacement in the X-direction yields 

K = -F x /e 0 , k = F Y /e 0 (8) 
Since this static displacement method has no dynamic motion, no damping 
or added-mass terms can be evaluated. 

A second approach to rotordynamic coefficient identification is 
utilized by Childs[15]. Depicted in Fig. 5, this method uses a circular 
orbit of the rotor within the seal. The rotor is mounted eccentrically 
on a shaft which rotates. Thus, the rotor precesses in a circular orbit 
at the same rate and direction as shaft rotation. This synchronous 
precession provides for the determination of the radial and tangential 
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components of the seal reaction force. The F r and F components are 
obtained through integration of the measured pressure distribution along - 
and around the seal periphery. Expressing measured rotor motion as 

X = e 0 cos(uit) 

(9) 

Y = e 0 sin(mt) 

for small circular orbit of radius e 0 and precessional frequency , 
and substituting into equation (2) yields the radial and tangential 
force coefficcient definitions 

F r / e 0 = Mm 2 - Cm ~ K = -K e f + M e f.ai 2 


( 10 ) 


F / e Q = k - Cu> = -C e f.w 

where the cross-coupled mass coefficient is assumed negligible with 
respect to the influence of k and C. Because the cross-coupled 
coefficients k and c are linear functions of u, identification of the 
individual dynamic coefficients is not possible in this method. However, 
equivalent direct stiffness, damping, and added-mass coefficients can be 
calculated as indicated in equation (10). 

Independent rotation and precession control, as shown in Fig. 6, is 
a third testing method which is currently employed both in impeller and 
seal studies [16], [17], [18]. Various means are used to produce a 
circular orbit (precession) of the rotor or impeller at a rate different 
from its rotational speed. For a small circular orbit of radius e Q and 
precessional frequency Q , the measured precessional motion of the rotor 
is 


X = e Q cos (fit) 

Y = e 0 sin(flt) 

The Fy and Fy reaction force components are measured and can be 




expressed 
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. F X = F X s sin(flt) + F xc cos(Qt) 

( 11 ) 

Fy = Fys sin(nt) + Fyc cos(Qt) 

By substituting these expressions into equation (2) and equating 
coefficients of sine and cosine terms, the following equations are 
obtained 


" F XC / e 0 = K + cn - Mn 2 
- F X s / e 0 = k - Cfl - mfl 2 


- Fyc / ©o = + Cfi + mft 2 

" F YS / e 0 = K + c.n + MQ 2 


( 12 ) 


Hence, by measuring the reaction force components and rotor motion at 
two different precession frequencies, eight equations in six unknowns 
are obtained, and the rotordynamic coefficients can be calculated. 

A fourth testing method has been used by lino and Kaneko [19] for 
determining dynamic coefficients, and this same method is employed at 
the TAMU gas seal test facility. An external hydraulic shaker is used 
to impart translatory harmonic motion to the rotating seal, and rotor 
motion relative to the stator and the reaction force components acting 
on the stator are measured. 


Fig. 7 shows the manner in which the rotor could be positioned and 
oscillated in order to identify the dynamic coefficients of the seal for 
small motion about e Q . If the added-mass terms are assumed negligible, 
equation (1) is rewritten 




C XX<Eo) C XY ( £ 0 ) 
C Yx(e 0 ) C YY(£q ) 



( 13 ) 


First, harmonic horizontal motion of the rotor is assumed, where 
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X = e 0 + A ain(nt) + B cos(Qt) 
X = Aft cos(nt) - BQ sin(nt) 

Y = Y = 0 


This yields small motion parallel to the static eccentricity vector, 
where n is the shaking frequency. In a similar fashion, the X and 
Y-direction force components can be expressed 

F X = F XS sin(Qt) + Fxc cos(flt) 


(14) 

Ff = Fys sin(nt) + Fyc cos(flt) 

Substituting these expressions into equation (13) and equating 
coefficients of sine and cosine terms yields the following four 
equations 


F XS= K XX A " C XX B 


F XC = K XX B + C XX A 
f YS = k YX A " C YX 8 
f YC = k YX 8 + C Y x A 


(15) 


Solving this system of four equations in four unknowns defines the 
dynamic coefficients as 

Kxx(^o) = ( F XC B + F XS A ) / (A 2 + B 2 ) 

Kyx(^o^ = ( F YS a + f YC B) / (A 2 + B 2 ) 

( 16 ) 

Cxx(^o) = ( F XC A - F XS B) / Q( A 2 + B 2 ) 

Cyx(£ 0 ) = ( F YC A - Fys B) / fl(A 2 + B 2 ) 


Therefore, by measuring the reaction forces due to known rotor 
motion, determining the Fourier coefficients (A, B, Fxs> fxc> F YS> F YC^> 
and substituting into the above definitions, the indicated dynamic 
coefficients can be identified. If the rotor is shaken about a centered 
position, then the process is complete. Since the linearized model has 
skew-symmetric stiffness and damping matrices, all of the coefficients 
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are identified. If, however, the rotor is shaken about an eccentric 
position as initially postulated, then it must be shaken vertically 
about that same point in order to complete the identification process. 

Assuming harmonic vertical motion of the rotor, as defined by 

X - e 0 , X = 0, 

Y = A sin(Qt) + B cos(flt), and 

Y = AQ cos(flt) - Bfl sin(Qt), 

yields oscillatory motion that is perpendicular to the assumed static 
eccentricity vector. A similar process as before results in the 
coefficient definitions 

Kyy(e 0 ) = < F XS A + F XC B ) / (A 2 + B 2 ) 

K XY( e o) = ~ ( f YC b + f YS A) / (A 2 + B 2 ) 

(17) 

C Y Y(e 0 ) = < f XC A - F XS B) / «(A 2 + B 2 ) 

Cxy( £ o) = ( f YS B - Fyc A) / n(A 2 + B 2 ) 

All eight dynamic coefficients are thus determined by alternately 
shaking the rotor at one frequency £2 in directions which are parallel 
and perpendicular to the static eccentricity vector. 



TEST APPARATUS OVERVIEW 


Detailed design of the TAMU gas seal apparatus was carried out by 
J.B. Dressman of the University of Louisville. It is of the external 
shaker configuration, and the dynamic coefficient identification process 
is as described in the latter part of the preceding section. 

Considering both the coefficient identification process and 
Nelson's analysis, some objectives for the design of the test apparatus 
are apparent. First, in order to determine the dynamic coefficients, 
the apparatus must provide for the necessary rotor motion within the 
seal, and measurement of the reaction-force components due to this 
motion must be possible. Secondly, it would be advantageous (for 
purposes of comparison) if the apparatus could provide the same variable 
seal parameters afforded by Nelson's analysis ( i . e . , pressures, seal 
geometry, rotor rotational speed, fluid prerotation, and rotor/stator 
surface roughness). With this capability, the influence of each 
independent parameter could be examined and compared for correlation 
between theoretical predictions and experimental results. 

With these design objectives in mind, the discussion of the test 
apparatus is presented in three sections. The first section. Test 
Hardware, describes how the various seal parameters are physically 
executed and controlled. For example, the manner in which the dynamic 
"shaking" motion of the seal rotor is achieved and controlled is 
described in this section. The second section, Instrumentation, 

describes how these controlled parameters, such as rotor motion, are 
measured. Finally, the Data Acquisition and Reduction section explains 
how these measurements are used to provide the desired information. 
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TEST HARDWARE 

This section deals only with the mechanical components and 
operation of the test apparatus. It is intended to provide answers to 
the following questions: 

1 ) How is the static position of the seal rotor controlled? 

2) How is the dynamic motion of the rotor executed and controlled? 

3) How is compressed air obtained and supplied to the apparatus, 

and how is the pressure ratio across the seal controlled? 

4) How is the incoming air prerotated before it enters the seal? 

5) How are the seal rotor and stator mounted and replaced? 

6) How is the seal rotor driven (rotated)? 

Recalling the rotordynamic coefficient identification process 
described earlier, the external shaker method requires that the seal 
rotor be set in some static position and then be oscillated about that 

point. The test apparatus meets those requirements by providing 

independent static and dynamic displacement control, which are described 
below. 

Static Displacement Control . The test apparatus is designed to provide 
control over the static eccentricity position both horizontally and 
vertically within the seal. The rotor shaft is suspended 
pendulum-fashion from an upper, rigidly mounted pivot shaft, as shown in 
Figs. 8 and 9. This arrangement allows a side-to-side (horizontal) 
motion of the rotor, and a cam within the pivot shaft allows vertical 
positioning of the rotor. 

The cam which controls the vertical position of the rotor is driven 
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Fig. 8 Components used for static and dynamic displacement of seal rotor. 
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Fig. 9 Test apparatus. 
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by a remotely-operated DC gearhead motor, allowing accurate positioning 
of the rotor during testing. Horizontal positioning of the rotor is 
accomplished by a Zonic hydraulic shaker head and master controller, 
which provide independent static and dynamic displacement or force 
control. The shaker head is mounted on an I-beam support structure, and 
can supply up to 4450 N (1000 lbf) static and 4450 N dynamic force. As 

illustrated in Fig. 8, the shaker head output shaft acts on the rotor 

shaft bearing housing, and works against a return spring mounted on the 
opposite side of the bearing housing. The return spring maintains 
contact between the shaker head shaft and the bearing housing, thereby 
preventing hammering of the shaker shaft and the resulting loss of 
control over the horizontal motion of the rotor.. 

Dynamic Displacement Control . The dynamic motion of the seal rotor 

within the stator is horizontal. In addition to controlling the static 
horizontal position of the rotor, the Zonic shaker head moves the rotor 
through horizontal harmonic oscillations as the test is run. A Wavetek 
function generator provides the sinusoidal input signal to the Zonic 
controller, and both the amplitude and frequency of the rotor 

oscillations are controlled. 

Although the test rig design provides for dynamic motion of the 
rotor only in the horizontal X-direction, all of the coefficients for 
either seal model (equation (3) or (13)) can still be determined. As 
Fig. 10 shows, the required rotor motion perpendicular to the static 
eccentricity vector can be accomplished in an equivalent manner by 
statically displacing it the same amount (e 0 ) in the vertical direction 
and continuing to shake horizontally. 


In addition to providing control over the rotor's static position 
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and dynamic motion, the test apparatus allows other seal parameters to 
he controlled independently, providing insight into the influence these 
parameters have on seal behavior. These parameters coincide with the 
variable input parameters for Nelson's analysis, and they include: 

1 ) pressure ratio across the seal , 

2) prerotation of the incoming fluid, 

3) seal configuration, and 

4) rotor rotational speed. 

Pressure Ratio . The inlet air pressure and attendant mass flow rate 
through the seal are controlled by an electric-over-pneumatically 
actuated Masoneilan Camflex II flow control valve located upstream of 
the test section. An Ingersoll-Rand SSR-2000 single stage screw 
compressor rated at 34 m 3 /min @ 929 kPa (1200 scfm 120 psig) provides 
compressed air, which is then filtered and dried before entering a surge 
tank. Losses through the dryers, filters, and piping result in an 
actual maximum inlet pressure to the test section of approximately 
722 kPa (90 psig) and a maximum flow rate of 27 mVmin (950 scfm). A 
four-inch inlet pipe from the surge tank supplies the test rig, and 
after passing through the seal, the air exhausts to atmosphere through a 
manifold with muffler. 

Inlet Circumferential Velocity Control . In order to determine the 
effects of fluid rotation on the rotordynamic coefficients, the test rig 
design also allows for prerotation of the incoming air as it enters the 
seal. This prerotation introduces a circumferential component to the 
air flow direction, and is accomplished by guide vanes which direct and 
accelerate the flow towards the annulus of the seal. The vanes 


are machined from brass disks, and Fig. 11 illustrates the vane 
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configuration. Three sets of guide vanes are available; one rotates the 
flow in the direction of rotor rotation, another introduces no fluid 
rotation, and the third rotates the flow opposite the direction of rotor 
rotation. 

Seal Configuration . The design of the test rig permits the installation 
of various rotor/stator combinations. As shown in Figs. 12-15, the 
stator is supported in the test section housing by three Kistler quartz 
load cells in a trihedral configuration. Figs. 12 and 13 show the 
smooth-rotor/smooth-stator seal, while the smooth-rotor/honeycomb-stator 
seal is illustrated in Figs. 14 and 15. The seal rotor is press-fitted 
and secured axially by a bolt circle to the rotor shaft. Seals with 
different geometries (i.e., clearances, tapers, lengths) can be tested, 
as well as seals with different surface roughnesses. The replacement of 
these rotor/stator combinations can be accomplished with minimal 
downtime. 

Rotational Speed . A Westinghouse 50-hp variable-speed electric motor 
drives the rotor shaft through a belt-driven jackshaft arrangement. 
This shaft is supported by two sets of Torrington hollow-roller bearings 
[20]. These bearings are extremely precise, radially preloaded, and 
have a predictable and repeatable radial stiffness. Axial thrust due to 
the pressure differential across the seal is absorbed by a flat, 
roller-type, caged thrust bearing at the rear of the rotor. Both the 
shaft and thrust bearings are lubricated by a positive-displacement 
gear-type oil pump. 

Different jackshaft drive-pulleys can be fitted to provide up to a 
4:1 speed increase from motor to rotor shaft, which would result in 
a rotor shaft speed range of 0-21,200 cpm. Current design limitations, 
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Fig. 12 Cross sectional view of test section showing smooth stator. 
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Fig. 14 Cross-sectional view of test section showing honeycomb stator. 
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however, prevent the attainment of this upper rotational speed. High 
bearing temperatures, reduction of interference in the rotor-rotor shaft 
fitment due to inertia-induced radial growth of the rotor inside 
diameter, and excessive stresses in the drive-pulleys have served to 
limit shaft speed. The highest rotational speed attained at the time of 
this writing is 8500 cpm, although design modifications to allow higher 
speeds are under investigation. 

To conclude this discussion of the test hardware, two views of the 
complete test apparatus are included. Fig. 16 shows the assembled rig, 


while an exploded view is provided in Fig. 17. 
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Fig. 16 Test apparatus assembly. 
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INSTRUMENTATION 

Having discussed what seal parameters can be varied, and how the 
variations are implemented, the measurement of their respective effects 
can now be described. The types of measurements which are made can be 
grouped into three categories: 

1 ) rotor motion, 

2) reaction-force measurements, and 

3) fluid flow measurements. 

These categories are described individually in the sections that follow. 
Rotor Motion Measurements . The position of the seal rotor within the 
stator is monitored by two Bently-Nevada eddy-current proximity probes, 
mounted in the test section housing. These probes are located 90 degrees 
apart, and correspond to the X and Y- directions. The proximity probes 
are used to determine the static position and dynamic motion of the 
rotor, and their resolution is 0.0025 mm (0.1 mil). 

Reaction-Force Measurements . Reaction forces arise due to the static 
position and dynamic motion of the seal rotor within the stator. The 
reaction forces (F)(, Fy) exerted on the stator are measured by the three 
Kistler quartz load cells which support the stator in the test section 
housing. When the rotor is shaken, vibration is transmitted to the test 
section housing, both through the thrust bearing and through the housing 
mounts. The acceleration of the housing and stator generates unwanted 
inertial "ma" forces which are sensed by the load cells, in addition to 
those pressure forces developed by the relative motion of the seal rotor 
and stator. For this reason, PCB piezoelectric accelerometers with 
integral amplifiers are mounted in the X and Y-directions on the stator, 
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as shown in Figs. 13 and 15. These accels allow a (stator mass) x 
(stator acceleration) subtraction to the forces (Fx, Fy) indicated by 
the load cells. With this correction, which is described more fully in 

the next section, only the pressure forces due to relative seal motion 

are measured. 

Force measurement resolution is a function of the stator mass and 
the resolution of the load cells and accelerometers. Accelerometer 
resolution is 0.005 g, which must be multiplied by the stator mass in 
order to obtain an equivalent force resolution. The masses of the 

stators used in the test program reported here are 11.4 kg(25.2 lb) and 
3-9^ kg(8.69 lb), corresponding to the smooth and honeycomb stators, 
respectively. Hence, force resolution for the accelerometers is 0.560 N 
(0.126 lb) and 0.191 N (0.043 lb), for each stator, respectively. 

Resolution of the load cells is 0.089 N(0.02 lb). Therefore, the 
resolution of the force measurement is limited by the accelerometers . 
With a stator with less mass, and/or accelerometers with greater 
sensitivity, force resolution could be improved. 

Fluid Flow Measurements . Fluid flow measurements include the leakage 
(mass flow rate) of air through the seal, the pressure gradient along 
the seal axis, the inlet fluid circumferential velocity, and the 
entrance pressure loss. 

Leakage is measured with a Fischer & Porter vortex flowmeter 
located in the piping upstream of the test section. Resolution of the 
flowmeter is 0.0014 m 3 (0.05 acf ) , and pressures and temperatures up and 
downstream of the meter are measured for mass flow rate determination. 

For measurement of the axial pressure gradient, the stator has 
pressure taps drilled along the length of the seal in the axial 
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direction. These pressures, as well as all others, are measured with 
a 0-1.034 MPa (0-150 psig) Scani valve differential-type pressure 
transducer through a 48 port, remotely-controlled Scani valve model J 
scanner. Transducer resolution is 0.552 kPa (0.08 psi). 

In order to determine the circumferential velocity of the air as 
it enters the seal, the static pressure at the guide vane exit is 
measured. This pressure, in conjunction with the measured flowrate and 
inlet air temperature, is used to calculate a guide vane exit Mach 
number. A compressible flow continuity equation 

m = P ex A ex M ex C(Y/R g T t ) (1 + (Y-1)M ex 2 / 2)] 1/2 (18) 
is rearranged to provide a quadratic equation for Mex 

M ex 2 = M + 1 ♦ 4 ( (Y-1 )/2Y ) (m R g T t / p ex A ex ) 2 } / (Y-1 ) (19) 
where Y is the ratio of specific heats and R g is the gas constant for 
air, T t is the stagnation temperature of the air, p ex is the static 
pressure at the vane exit, and A ex is the total exit area of the guide 
vanes. Since all of the variables in the equation are either known or 
measured, the vane exit Mach number, and therefore the velocity, can be 
found. 

In order to determine the circumferential component of this inlet 
velocity, a flow turning angle correction, in accordance with Cohen[21], 
is employed. The correction has been developed from guide vane cascade 
tests, and accounts for the fact that the fluid generally is not turned 
through the full angle provided by the shape of the guide vanes. With 
this flow deviation angle calculation, the actual flow direction of the 
air leaving the vanes (and entering the seal) can be determined. Hence, 
the magnitude and direction of the inlet velocity is known, and the 
appropriate component is the measured inlet circumferential velocity. 
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The entrance pressure-loss coefficient, defined in equation (6), is 
determined from the measured pressures just upstream of and just inside 
the seal. An entrance Mach number is calculated in the same manner as 
outlined previously, using the measured pressure immediately inside the 
seal and the annular area between the rotor and stator. This entrance 
Mach number, and the ratio of the seal entrance/ guide vane exit 
pressures are substituted into equation (6), and the entrance loss 
coefficient, k, is determined. 
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DATA ACQUISITION AND REDUCTION 

With the preceding explanations of how the seal parameters are 
varied, and how these parameters are measured, the discussion of how the 
raw data is processed and implemented can begin. Data acquisition is 
directed from a Hewlett-Packard 9816 ( 1 6 — bi t ) computer with disk drive 
and 9.8 megabyte hard disk. The computer controls an H-P 6940B 
multi programmer which has 12-bit A/D and D/A converter boards and 
transfers control commands to and test data from the instrumentation. 

As was previously stated, the major data groups are seal 
motion/reaction force data and fluid flow data. The motion/reaction 
force data are used for dynamic coefficient identification. The 

hardware involved includes the load cells, accelerometers', X-direction 
motion probe, a Sensotec analog filter unit, a tunable bandpass filter, 
and the A/D converter. The operation of these components is illustrated 
in Fig. 18, and their outputs are used in a serial sampling scheme which 
provides the computer with the desired data for reduction. Recalling 
the discussion of the reaction force measurements in the preceding 
section, a (stator mass) x (stator acceleration) subtraction 
from the Indicated load cell forces is necessitated due to vibration of 
the stator and test section housing. This subtraction is performed 
with an analog circuit, and results in corrected Fy and Fy force 
components due to relative seal motion. 

The forced oscillatory shaking motion of the seal rotor is the key 
to the operation of the serial synchronous sampling (SSS) routine which 
is employed. The frequency of the rotor oscillation is set by a function 
generator, and rotor motion is sensed by the X-direction motion probe. 
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Fig. 18 Signal conditioning schematic for data acquisition. 
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The motion signal is filtered by the narrow bandpass filter, and is used 
as a trigger signal for the SSS routine. Upon the operator's command, 
the SSS routine is enabled, and the next positive-to- negative crossing 
of the filtered motion signal triggers a quartz crystal clock/timer. 
Ten cycles of the corrected Fx(t) signal are sampled, at a rate of 100 
samples/cycle . The second positi ve-to-negative crossing of the filtered 
motion signal triggers the timer and initiates the sampling of ten 
cycles of the Fx<t) signal. Finally, the third positive-to-negati ve 
crossing triggers the timer again, and ten cycles of the corrected X ( t ) 
signal are sampled. Thus, at every test condition, 1000 data points are 
obtained for Fx(t^ ) ,Fy(ti ) , and X(t^), and the data arrays are stored in 
computer memory. 

Some important points need to be stressed concerning this 
force/motion data acquisition. First, the bandpass filter is used only 
to provide a steady signal to trigger the timer/clock. Any modulation 
of the motion signal due to rotor runout is eliminated by this filter, 
as long as the rotational frequency and shaking frequency do not 
coincide. Therefore, the shaking frequencies are selected to avoid 
coincidence with running speeds. However, the rotor motion and corrected 
force signals which are sampled and captured for coefficient 
identification are filtered only by a low-pass filter (500 Hz cutoff), 
and the effects of runout as well as shaking motion are present in the 
recorded data. A second point worth noting is that the sample rate is 
directly dependent on the shaking frequency. As the shaking frequency 
is increased, the sample rate (saraples/second ) also increases. In order 
to get the desired 100 samples/ cycle, shaking frequencies must be chosen 


to correspond to discrete sample rates which are available. 


Hence, the 
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frequency at which the rotor is shaken is carefully chosen to provide 
the desired sampling rate and a steady trigger signal. 

Most of the fluid flow data are used for the input parameters 
required by Nelson’s analysis. The upstream (reservoir) pressure and 
temperature, downstream (sump) pressure, and the inlet circumferential 
velocity (determined as outlined earlier) are provided directly. The 
friction-factor values of the rotor and stator are supplied in the 
form of Hirs coefficients, which are obtained from the pressure 
distribution data in the manner described below. 

Recalling the discussion of Hirs' turbulent bulk-flow fluid model, 
the model assumes that the wall shear stresses can be written as in 
equation (4). For the gas seals discussed here, an adiabatic, 
compressible flow with friction analysis is employed, and the measured 
pressure gradient and mass flow rate (leakage) data are used to 
calculate a friction factor coefficient. A, for each test condition. 
From the A versus R a and w data, the Hirs coefficients mr , nr of the 
friction factor formula 

A = nr R a mr [ 1 + 1 / 4b 2 ] ( 1 + n>r)/2 f b - U / Ru (20) 

are calculated on a least-square basis. For the 
smooth-rotor/smooth-stator combination, the values are assumed to apply 
for both the rotor and stator. Hence, for this case, mr=ms and nr=ns. 

For the smooth-rotor/honeycomb-stator combination, a combined A is 
measured, which is related to the rotor A r and (honeycomb) stator A s by 

A c = (A r + A 3 ) / 2 (21 ) 

and hence, 


( 22 ) 
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Therefore, A s is determined from measured data for A c and a 
calculated value for A r from Equation (21) with experimentally- 

determined values for mr and nr. Then, as before, the A s versus R a and w 
data are used to calculate the Hirs coefficients for the honeycomb 
stator . 

As stated previously, the Hirs coefficients for the seal rotor and 
stator are required input parameters for Nelson's analysis, as are the 
fluid flow conditions up and downstream of the seal and the rotational 
speed of the rotor. The appropriate input parameters for each specific 
test case can be provided for Nelson's analysis from static test results 
and measurements. In this manner, a point-by-point comparison of 
theoretical predictions to experimental results can be made for leakage 
through the seal, axial - pressure distribution, entrance - loss 
coefficient, and rotordynamic coefficients. 
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TEST PROCEDURES 

At the start of each day's testing, the force, pressure, and 
flowmeter systems are calibrated. The total system, from transducer to 
computer, is calibrated for each of these variables. The force system 
calibration utilizes a system of pulleys and known weights applied in 
the X and Y-directions . An air-operated dead-weight pressure tester is 
used for pressure system calibration, and flowmeter system calibration 
is achieved with an internal precision quartz clock which simulates a 
known flowrate. 

All of the tests performed to date have been made with the rotor 
executing small motion about a centered position. A typical test begins 
by centering the seal rotor in the stator with the Zonic hydraulic 
shaker, starting airflow through the seal, setting the rotational 
speed of the rotor, and then beginning the shaking motion of the rotor. 
Data points are taken at rotational speeds of 200, 500, and 1000-8000 
cpm, in 1000 cpm increments. At each rotational speed, the inlet 
pressure is varied and data points are taken at one unchoked flow and 
four choked flow conditions. For each test case (i.e., one particular 
running speed, shaking frequency, inlet pressure, and prerotation 
condition), the measured leakage, rotordynamic coefficients, axial 
pressure distribution, and entrance loss coefficient are determined and 
recorded. 

This test sequence is followed for each of three different shaking 
frequencies, and for three inlet swirl directions (with rotor rotation, 
opposite rotation, and no rotation). Therefore, fifty data points are 
taken per test (i.e. one shaking frequency and inlet swirl combination), 
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with a total of nine tests (for small motion about a centered 
position) made per seal. 
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RESULTS 

The test results reported here were developed as part of an 
extended, joint NASA-USAF funded research program for annular gas seal 
studies. Tests were carried out on a smooth-rotor/smooth-stator seal 
and a smooth-rotor/honeycomb-stator seal. The dimensions and pertinent 
data for each are included in Table 1. The honeycomb stator insert, 
provided by the Rocketdyne division of Rockwell International, is the 
turbine interstage seal of the HPOTP (High-Pressure Oxygen Turbopump) of 
the SSME (Space Shuttle Main Engine). Fig. 19 illustrates the honeycomb 
configuration. 

The test program had the following objectives: 

1) Acquire leakage, friction factor, and entrance-loss data for 
smooth and honeycomb seals. 

2) Compare predictions from current theory to test results. 

3) Compare the stability performances of a smooth-rotor/ 
honeycomb-stator and a smooth constant-clearance annular seal. 

The Hirs coefficients for both seals were determined in the manner 
described previously. The values of these coefficients are listed in 
Table 2. Relative roughness values based on measured Hirs coefficients 
as suggested by Colebrook [22] are also included. Colebrook's 
formulation, 

(4 nsR a ms r 1/2 = -2 log ( ( ( e/2C ) / 3-7) + (2.51 / (R^4 nsRams ))) (23) 

is a curve-fit of experimental data obtained for fluid flow through 
pipes with various wall roughnesses. The appropriate stator 

coefficients are substituted to obtain the relative roughness (e/2C) 
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Table 1. Test seal specifications. 

Smooth-rotor/ Smooth-rotor/ 

smooth stator honeycomb stator 

Rotor 

Diameter: 15.136 cm (5.959 in) 14.453 cm (5.690 in) 

Material: 304 Stainless steel 304 Stainless steel 

Surface 

roughness: 0.102 pm (4 pin) 0.127 um (5 yin) 

Stator 

Diameter: 15.283 cm (6.017 in) 14.614 cm (5.754 in) 

Material: 304 Stainless steel 6061-T6 Aluminum 

Surface 

roughness: 0.140 ym (5.5 yin) 1.575 mm (0.062 in) 

Comb 

Radial clearance: 0.7366 mm (29 mil) 0.8065 mm (31-8 mil) 


Seal length: 


5.080 cm (2.00 in) 


2.540 cm (1 .00 in) 



Fig. 19 Honeycomb stator insert detail 
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Table 2. Friction-factor data. 

Smooth-rotor/ 
smooth-stator 

Hirs coefficien ts 

Rotor ns: 

Rotor ms: 

Stator ns: 

Stator ms: 


1.44x10" 4 1.44x10 -4 

1.44x10 _i< 4.93xl0 _JJ 

(Note: The relative roughness values shown are averages over an axial 

Reynolds number range of 20,000 - 80,000.) 



0.187 

0.187 

-0.333 

-0.333 

0.187 

0.187 

-0.333 

-0.0778 


Smooth-rotor/ 

honeycomb-stator 
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values. It should be noted that friction factor data for honeycomb 
seals have not been previously published. 

The results provided here are grouped in static (leakage, pressure 
distribution, entrance loss coefficient) and dynamic (rotordynamic 
coefficient) sections. A one-to-one comparison of the smooth and 
honeycomb seal configuration is precluded, however, due to differences 
in seal length, nominal clearance, and inlet guide vane configuration, 
as illustrated in Figs. 12 and 14. 

Static Results . Figs. 20 and 21 and Table 3 show a comparison of the 
theoretical and experimental leakage through the seal for various fluid 
prerotation conditions. The figures show the leakage at various 
pressure ratios (reservoir pressure / sump pressure). In the table, the 
leakage has been averaged over all speeds and pressure ratios, and is 
presented in ratio form (Theory/Exper iment ) . The comparison shows that 
for both the smooth and the honeycomb seal, leakage is underpredicted 
for the non-prerotated case. Conversely, for fluid prerotation either 
in or opposing the direction of rotor rotation, the leakage is 
overpredicted for both seals. The maximum error is approximately 7 . 5 %, 
occurring for the smooth seal with prerotation in the direction of rotor 
rotation. Average error for the smooth seal is 1.7%, and for the 
honeycomb seal is \%. 

The pressure gradient plots are included in Figs. 22-28. Fig. 22 
illustrates the negligible effect of running speed on the pressure 
distribution in the seal. This plot has ten curves (corresponding to 
the ten rotational speed increments) plotted. This accounts for the 
heavy lines which appear in some cases. This particular plot is of 
the experimental data for the non-prerotated smooth seal case, however, 
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VERSUS PRESSURE RATIO 



Fig. 21 Honeycomb seal leakag 
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Table 3. Theory versus experiment leakage comparison. 
(Theory/Experiment ) 


Fluid prerotation 
direction 


Smooth-rotor/ Smooth-rotor/ 

smooth-stator honeycomb- stator 


With rotor rotation 1.075 1.047 

(standard deviation) 0.012 0.008 

No prerotation 0.9684 0.9712 

(standard deviation) 0.012 0.005 

Opposite rotor 

rotation 1.007 1.013 

(standard deviation) 0.012 


0.004 



FILES OF TYPEc SQ1D11C 7 
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Fig. 22 Axial pressure gradient in smooth seal with no prerotation, unchoked flow 
ten running speeds. 
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none of the pressure plots show any appreciable variation due to running 
speed. 

Figs. 23-28 show the theoretical and experimental pressure data for 
each of the seals under various prerotation conditions. Due to the 
absence of running speed dependence, only one speed is plotted for each 
inlet pressure condition. The numbers on the plotted lines refer to the 
inlet pressure, where 1 corresponds to approximately 186.9 kPa (12.4 
psig), and 2 through 5 correspond to 308.2 kPa (30 psig), *446. 1 kPa (50 
psig), 584.0 kPa (70 psig), and 721.9 kPa (90 psig), respectively. The 
lowest pressure corresponds to unchoked flow through the seal, while the 
others are choked. The shapes of the pressure-gradient plots show 
fairly good correspondence between theory and experiment. This to be 
expected, however, since the Hirs' coefficients used in the analysis 
come directly from the measured pressures. The best agreement occurs 
for the non-prerotated flow in both the smooth and honeycomb seals. For 
prerotated flow in either direction, the theoretical gradient is shifted 
up slightly for both seal configurations. This upward shift is partly 
due to a total pressure correction that is made. When the flow is 
prerotated by the guide vanes, it is accelerated as well as turned, and 
the measured static pressure at the vane exit decreases. This explains 
why the experimental plots show lower seal entrance pressures for either 
prerotated case than for the non-prerotated case. Nelson's analysis, 
however, assumes that the supply pressure upstream of the seal is the 
total pressure. Hence, the axial component of the fluid velocity as it 
leaves the guide vanes is used to calculate an effective total pressure, 
which is higher than the measured static pressure. This corrected 
pressure is then input as the reservoir pressure to Nelson's analysis. 
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Fig. 23 Axial pressure gradient in smooth seal with prerotation in the direction of rotor rotation 
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The entrance loss coefficient, k, also may have some bearing on the 
upward shift seen in some of the pressure gradient plots. Plots of 
(k+1) versus axial Reynolds number are included in Figs. 29 — 3 ^ • Ten 

experimental curves, corresponding to running speed are plotted on each. 
Recalling Deissler's curve fit employed by Nelson and plotted in Fig. 3, 
experimental results show loss terms (k+1 ) significantly higher than 
those predicted for both prerotated honeycomb seal cases. Agreement 
between theory and experiment is fairly good for the non- prerotated 
cases for both seals. For the ■ smooth seal with prerotation in the 
direction of rotor rotation, the loss coefficient is overpredicted, with 
the experimental results indicating a negative k. 

Dynamic Results . Dynamic tests were performed at shaking frequencies of 
58.8, 74.6, and 124.6 Hz. As was discussed in the Data Acquisition 
section of this report, these frequencies were chosen to provide the 
desired sample rate and a steady trigger signal. The dynamic 
coefficients obtained at the two lower frequencies are essentially the 
same. At the 124.6 Hz shaking frequency, however, correspondence of the 
data to that obtained at the lower frequencies is unsatisfactory. In 
seeking to explain the discrepancy, tests were run to determine the 
relative transfer function of the test apparatus. The plots in Fig. 35 
show the results of these tests, and indicate a resonance of 

the apparatus occurring at approximately 25 Hz (the drop in phase 

difference at approximately 45 Hz corresponds to a resonance of the 
shaker support structure). As the shaking frequency is increased above 
this, the input force levels required to achieve a given motion 
amplitude increase rapidly. At the 124.6 Hz shaking frequency, 

attainable motion amplitude is about 50% of that achieved at the 
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Fig. 33 Entrance-loss for honeycomb seal 
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Fig. 35 Relative transfer function of test apparatus 
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58.8 and 74.6 Hz frequencies. Therefore, one possible explanation for 
the poor agreement between the results is that as motion amplitude 
decreases, so does the force measured by the load cells, and the 
measurement system resolution suffers. 

The plots of the rotordynamic coefficients are found in Figs. 
36-43. These plots include both the theoretical and experimental data. 
The coefficients are plotted versus the reservoir / sump pressure ratio, 
and the solid lines correspond to the theoretical data. The 
experimental results are indicated by symbols. The location of the 
symbols represents the average value of the coefficient (averaged 
over all of the running speeds) at each particular inlet pressure 
condition, and the vertical lines drawn through the symbols signify the 
range over which they varied through the speed range. The test results 
plotted here were obtained by shaking the rotor with an amplitude of 
seven mils at 74.6 Hz. 

Dynamic Results - Smooth Seal . For the smooth seal, direct stiffness 
(Fig. 36) is overpredicted for the non-prerotated condition, and 
underpredicted for both prerotated conditions. Best agreement is seen 
in the case for prerotation opposite rotor rotation, and the trend of 
increasing stiffness with increasing pressure ratio compares favorably. 

In the cross-coupled stiffness comparison for the smooth seal 
(Fig. 37), theory overpredicts for both prerotation conditions, and 
underpredicts for the straight flow case. In every instance, agreement 
becomes worse with increasing pressure ratio. The non-prerotated case 
shows a divergence both in magnitude and sign. It should be noted, 
however, that the magnitudes for this case are significantly smaller 


than for either of the prerotated cases. 







SMOOTH SEAL EXPERIMENT & THEORY 
VERSUS PRESSURE RATIO 





Agreement between theory and experiment for direct damping (Fig. 38 ) 
is the most favorable of all the dynamic coefficients. Direct damping 
for the smooth seal shows an increase for increasing pressure ratio 
across the seal, with theory overpredicting slightly for the case of 
prerotation in the direction of rotor rotation. For the other 

prerotation conditions, the direct damping is slightly overpredicted at 
higher pressure ratios. 

Cross-coupled damping (Fig. 39) for the smooth seal generally shows 
agreement in the trends for the theoretical and experimental results. 
For prerotation in and opposing the direction of rotor rotation, the 
theory underpredicts cross-coupled' damping magnitude by approximately 
503», but shows a sign consistent with the test data. For the 
non-prerotated case, the theory predicts coefficients so small as to be 
considered negligible. This is not inconsistent with the test results, 
however, as the magnitudes for this case are significantly smaller than 
for either prerotated case. 

Dynamic Results - Honeycomb Seal . The honeycomb seal data, in 
general, shows the same correspondence between theory and experiment as 
the smooth seal. A notable exception, however, is in the direct 
stiffness coefficient (Fig. 40). For both prerotated cases, the 
predicted stiffness decreases with increasing pressure ratio, while the 
measured stiffnesses increase. This same predicted decreasing trend is 
shown for the non-prerotated case at the lower pressure ratios. 

In the cross-coupled stiffness comparison (Fig.4l), the theory 
underpredicts the magnitudes, but correctly predicts the signs of the 
coefficients. For the non-prerotated case, the predicted stiffnesses 
are essentially zero. The relative magnitudes of the experimental 
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Fig. 38 Direct damping of smooth seal. 
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Fig. 39 Cross-coupled damping of smooth seal 
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Fig. 40 Direct stiffness of honeycomb seal 
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Fig. 41 Cross-coupled stiffness of honeycomb seal 
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results for this same case in comparison to either prerotated case are 
also quite small, however. For prerotation in the direction of rotor 
rotation, theoretical cross-coupled stiffnesses are approximately 25$ 
less than experimental ones. For counter prerotation, theory 
underpredicts by about 50$. 

With the exception of the non- prerotated case, agreement between 
theory and experiment is fairly good for the direct damping coefficients 
(Fig. 42) of the honeycomb seal. In the non-prerotated case, theory 

underpredicts< the coefficients by approximately 46$. The prerotated 
cases show agreement to within approximately 10$. 

Theoretical results for the cross-coupled damping coefficients 
(Fig. 43) of the honeycomb seal are small enough to be considered 

negligible. In every case, the theory underpredicts the coefficients by 
a wide margin. However, the trend of increasing magnitude with 

increasing pressure ratio, as well as the signs of the coefficients, 

agree . 

One method in which the dynamic coefficients of the smooth and 
honeycomb seals can be directly compared is through their respective 
non-dimensional whirl frequency ratios. Whirl frequency ratio is defined 

Whirl frequency ratio = k / Cft , 

where ft is the shaking frequency. This ratio is a measure of the 
destabilizing influence of the cross-coupled stiffness with respect to 
the stabilizing influence of direct damping. Plots of whirl 
frequency ratio versus running speed are included in Fig. 44. The smooth 
seal plot shows a small, positive whirl frequency ratio over most of 
the running speed range. The honeycomb seal plot, however, shows a 
negative whirl frequency ratio. The negative sign arises due to a 
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Fig. 42 Direct damping of honeycomb seal 
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Fig. 43 Cross-coupled damping of honeycomb seal. 
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negative cross-coupled stiffness. This negative k exerts a stabilizing 
influence, resulting in a force which acts in the same direction as the 


direct damping force. 



CONCLUSIONS 


A seal-test facility has been developed for the study of various 
types of gas seals. A method of determining rotordynamic coefficients 
has been established, and consistent, repeatable results have been 
obtained. After some initial failures in the test apparatus, 
reliability has been satisfactory, and a complete set of experimental 
results can be acquired in a matter of days. 

The experimental and theoretical results of the preceding section 
support the following conclusions: 

(a) Theoretical results for leakage are consistent with test 
results. Slightly higher leakage occurs for cases where the flow is not 
prerotated. Agreement between theory and experiment is satisfactory, 
with a maximum error of approximately 7 . 5 %. 

(b) Experimental and theoretical results for the pressure 
distributions and entrance-loss coefficients are relatively insensitive 
to running speed for the ranges (0-8500 cpm) and seals tested to date. 

(c) The entrance-loss relationship (Deissler [13]) employed by 
Nelson is inconsistent with the test results for cases where the fluid 
is prerotated. For all such cases except one, the entrance-loss 
coefficient is underpredicted. 

(d) In the test results for the honeycomb seal, the steep entrance 
pressure-loss seems to extend partially inside the seal. Also, the 
measured pressure at the exit of the seal generally equals the back 
pressure, rather than being greater, as is predicted by compressible 
flow theory for choked flow. These same phenomena do not occur for the 
smooth seal, implying that perhaps the effective length of the honeycomb 
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seal is less than its actual physical length. 

(e) Test results for the direct stiffness of both the smooth and 
honeycomb seals show much greater sensitivity to fluid prerotation than 
predicted by theory. Prerotation of the fluid (in either direction) 
results in measured direct stiffnesses which are significantly larger 
than for no prerotation. Theory predicts the direct stiffness to be 
relatively insensitive to fluid prerotation. 

(f) Theoretical predictions of the influence of fluid prerotation 
on cross-coupled stiffness and damping are consistent with the test 
results. In general, theory underpredicts the magnitudes of these 
cross-coupled coefficients, while correctly predicting their trends 
with respect to prerotation. 

(g) Agreement between theory and test results for the direct 
damping coefficients is favorable. For both the smooth and honeycomb 
seal, direct damping is largest for no fluid prerotation. 

(h) Over the speed range tested, none of the rotordynamic 
coefficients show appreciable sensitivity to the rotational speed of the 
rotor. This may be due to the lack of development of significant shear 
forces in the seal. It appears that running speeds above those attained 
to date may be necessary to produce significant shear force effects. 

(i) For the non-prerotated case, the smooth seal has a positive 
cross-coupled stiffness, while k for the honeycomb seal is negative. 
This negative cross-coupled stiffness, and hence negative whirl ratio, 
indicates that the stability performance of the honeycomb seal is 


more favorable than that of the smooth seal. 
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